Effect of random fluctuations on quantum spin-glass transitions at zero temperature 
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We study the effects of random fluctuations on quantum phase transitions by the energy gap 
analysis. For the infinite-ranged spin-glass models with a transverse field, we find that a strong 
sample-to-sample fluctuation effect leads to broad distributions of the energy gap. As a result, the 
linear, spin-glass, and nonlinear susceptibilities behave differently from each other. The power-law 
tail of the distribution implies a quantum Griffiths-like effect that could be observed in various 
random quantum systems. We also discuss the mechanisms of the phase transition in terms of 
the energy gap by comparing the Sherrington-Kirkpatrick model and random energy model, which 
demonstrate the difference between the continuous and discontinuous phase transitions. 
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The quantum phase transition that could occur by 
changing parameters in the Hamiltonian is characterized 
by the closing of the energy gap between the ground and 
first excited states. The emergence of the gap is due to a 
quantum effect and we can find various new types of tran- 
sitions that are absent in classical thermodynamic sys- 
tems [l[ . The situation becomes more complicated when 
we treat disordered systems. We need to theoretically 
take an average over realizations of disordered parame- 
ters. We expect that the self-averaging property holds in 
the thermodynamic limit and that the average model de- 
scribes a real system represented by a specific realization. 
However, at low temperatures, the sample-to-samplc fluc- 
tuations become important and lead to phenomena that 
can never be seen in clean systems. Actually, the gap 
vanishing point strongly depends on the random sample 
and we can hardly identify the phase transition point. 
In a disordered phase of the system, it is also known 
that physical quantities are affected by the presence of 
finite clusters of ordered state. This Griffiths singularity 
0, Q is widely known in random spin systems as a gen- 
eral mechanism that induces random fluctuation effects. 
A similar effect is known in disordered electron systems 
as anomalously localized states in delocalized phase 
which shows the universality of the phenomenon for dis- 
ordered systems. 

The Griffiths singularity in the transverse-field Ising 
spin-glass model was found numerically by Rieger and 
Young Q , and Guo et al Q . They found that the nonlin- 
ear susceptibility at zero temperature is divergent below 
a point within the quantum paramagnetic phase in two- 
and three-dimensional systems. It was argued that mag- 
netically ordered clusters in a sample gives a small gap, 
which makes the susceptibility large. The same effect 
was also found by Fisher [7j for the same model in one 
dimension. Although this effect was called the quantum 
Griffiths singularity, the formation of the ordered cluster 
is essentially the same as that in a classical case. It is 
expected to disappear at large dimensions since clusters 
easily form at low dimensions. 



It is well known that quantum systems can be formu- 
lated using path integral formalism. Then, the quantum 
effect can be represented by fluctuations in an extra di- 
mension of imaginary time and the averaging over disor- 
der induces correlations between different times. This is 
a common feature that can be seen in all quantum sys- 
tems with disorder. In this letter, using simple quantum 
spin-glass models, we consider a Griffiths-like mechanism 
that can persist even in infinite-ranged models without 
the Griffiths singularity. 

We mainly discuss the Sherrington-Kirkpatrick (SK) 
model [1J in a transverse field as one of the simplest quan- 
tum spin-glass model. The Hamiltonian is written as 
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where af' z are Pauli matrices at site i, Ja are random 



interactions with an average [Jj. 



and a variance 



[Jfj] = J /N, N is the site number, and T is the trans- 
verse field. Several analyses showed that this model at 
zero temperature has a second-order phase transition be- 
tween the spin-glass and quantum paramagnetic phases 
at T/J - 1.4 - 1.5 0-E1. The SK model is the two- 
body interacting one, and one can generalize it into p- 
body ones. The model at p — > oo is known to be 
equivalent to the random energy model (REM) Its 
quantum version including the transverse field can be 
solved exactly and shows a discontinuous transition at 
T/J = \/h72 - 0.83 [H. We treat both the SK model 
([T]) and the REM since it is interesting to know the na- 
ture of the phase transition from the energy spectrum. 
We numerically diagonalize their Hamiltonian matrices 
by the Lanczos method. The number of spins, N, is 
taken up to 16. The corresponding size of the matrix is 
given by 2 N . The ensemble average is taken over more 
than 20000 samples. 

In Fig. [U we show the results of the average energy 
gap. In the SK model, the gap vanishing point is es- 
timated as T/J ~ 0.6, which is much smaller than the 
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expected transition point. The gap vanishing point can 
be roughly estimated using the perturbation theory as 
follows. When J = 0, energy level is exactly given by 
E n = -(N - 2n)T with n = 0, 1, • • • , N and each level 
has the degeneracy Nl/(N — n)\n\. When J is included, 
these degeneracies are lifted to form broad distributions 
of energy levels. The first excited state belongs to the 
sector n = 1. In this sector, the matrix element of the 



N x N Hamiltonian is given by H t 



Ci) 



-(N-2)T-J iL 



The second term corresponds to the Hamiltonian of the 
Gaussian random matrix theory, and the energy levels 
are known to form a semicircle distribution [l6j . The 
width of the distribution is given by 4 J and the first ex- 
cited energy level is given by E\ ~ — (N — 2)T — 2J. 
The unperturbed ground state with the sector n = is 
given by E ~ — AT and the energy gap is obtained as 
A ~ 2(T — J). Thus, the degenerate point is roughly 
estimated as T = J which is modified by correlations be- 
tween different sectors. Up to the second order in the 
perturbation theory, the degenerate point is estimated 
as T/J ~ 0.85. Since the perturbation becomes worse at 
small transverse fields, we cannot effectively predict the 
degenerate point. However, it can be definitely said that 
the perturbative correction decreases the value of the de- 
generate point and the point never approaches the phase 
transition point at T > J. 

On the other hand, for the REM, we find a change in 
behavior around the phase transition point. The average 
gap becomes minimum around the point and the fluctu- 
ation also changes its behavior. This result is consistent 



with the analysis discussed in Ref. 1171 where the transi- 
tion point is supposed to be equivalent to the point where 
the average gap becomes zero. 
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FIG. 1: Average energy gap A = [Ei — Eo] between the 
ground and first excited states for the transverse SK model 
and REM (N — 16). The variance of the energy gap 
±y/[(Ei — -Bo) 2 ] — A 2 is shown as the error bar. 



The gap fluctuation in Fig. [TJ can be clearly seen from 
the gap distribution function in Fig. [2] For the SK model, 
we find a single peak distribution with a long tail, which 
implies an important role of the fluctuation. In the case 
of the REM, we obtain more complicated distributions 
that are dependent on T. For a large T, we observe a 
single peak on the right-hand side of the distribution. 




FIG. 2: Gap distribution function (N = 16). 



When r decreases, this peak decreases and another peak 
at the left-hand side appears. The difference between 
the heights of the peaks becomes smaller and changes its 
sign around the transition point. This behavior together 
with the result of the average gap implies a structural 
change between two different configurations at the tran- 
sition point. 

In the classical spin-glass theory, it is known that three 
types of susceptibilities, namely, linear yi . spin-glass Xsg, 
and nonlinear Xni> play important roles [la. Il9j. When 
r = 0, the linear susceptibility is represented by the spin- 
glass order parameter q = [(erf) 2 ] as xi — /3(1 — q), and 
the nonlinear susceptibility can be expressed by the spin- 
glass one as Xni = /3(Xsg — 2/3 2 /3). The phase transition 
characterized by the divergence of Xsg can be found by 
observing the behavior of Xni- 

These relations are changed in quantum systems owing 
to the fluctuation effect. Using the spectral representa- 
tion, we can write the linear susceptibility as the sum 



(q) 



The first term 



of two contributions: xi = Xi ' + Xi 
x[ c ^ — P{x — <z) corresponds to the classical part and can 
be treated by the static approximation [2(|. There is no 
classical counterpart for x[ q ' '. At zero temperature, we 
have 



(q) 
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where \n) is the eigenstate with the energy E n , |0) is the 



nondegenerate ground state, and <r z — 53i=i cr f- I n the 
same way, the quantum parts of x sg and Xni are 
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The singularity of the susceptibility comes from the quan- 
tum part. It occurs when the energy gap of the first ex- 
cited level A reaches zero. Generally speaking, Xh Xsg, 
and Xni depend on 1/A 1 ' 2 ' 3 , respectively, and the energy 
gap distribution P(A) determines their critical behavior. 
When A is small, P(A) is assumed to have the power- 
law form P(A) ~ A fe . Then, xi> Xsg, and Xni diverge at 
k < 0, 1, and 2, respectively. It should be stressed here 
that the point where Xni diverges could be different from 
the point where Xsg does. 




FIG. 3: Inverse of the average susceptibility (quantum part) 
for the SK model with N = 14. Inset: The exponent k of the 
energy gap distribution function (N = 8,10,12,14, and 16 from 
above to below). The gap exponent is fitted by -P(A) ~ A fc 
near the origin. 



We plot the numerical result of the susceptibilities and 
the exponent k in the gap distribution function for the 
SK model in Fig. [3l The result implies that the three sus- 
ceptibilities diverge at different points, which are roughly 
the same as the estimates from the exponent. Although 
the second term in Eq. (fj| is neglected in the calcula- 
tion of Xni, we find in smaller systems that the second 
term does not change the divergent point. We also find 
JX\ ~ 2 at a point where Xsg diverges, which is consistent 
with results of the analysis discussed in Ref. Il2t 

In order to make sure that the gap between the ground 
and first excited states determines the divergence of the 
susceptibilities, we must examine the effect of higher- 
order excited states. For that purpose, we study the 
Fourier representations of a real-time correlation func- 
tion including contributions from all states: 
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From this spectral function, we see for the SK model 
that the main contribution comes from TV levels of the 
sector n = 1 if F is sufficiently large. This function was 
analytically calculated perturbatively in ref. [21| to yield 
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FIG. 4: Local correlation function xi^) f° r the SK model 
(N — 16). The shaded areas are contributions of the first 
excited state. 



a semicircle form. When T decreases, a sharp peak near 
the origin grows and comes into contact with the vertical 
axis [11]. As shown in Fig. : T this peak is made from 

the first excited state only. We have the relation X\ = 
/ dwx(w)/w, and the behavior of x( w ) around the origin 



determines the divergence of x- Actually, the 



divergence of the quantity J dujxiu) / ui 3 was used to find 
a phase transition [12j. Thus, we conclude that the first 
excited state determines the critical behavior. 
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FIG. 5: Inverse participation ratio (N=12). The summation 
over n in Eq. ([6]) is restricted to n = 0(1) for a plot with 
the caption n — 0(1). The error bars are smaller than the 
symbols. 



The above analysis cannot be applied to the REM. In 
this model, xi an d Xni do not diverge. They discontinu- 
ously jump at the transition point and stay constant in 
the spin-glass phase. The difference can be understood 
from the analysis of the inverse participation ratio de- 
fined as 



(6) 



This quantity is estimated as the inverse of the number of 
excited states included in of |0). A similar quantity was 
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used in disordered electron systems to study the localiza- 
tion 22|. As shown in Fig.[5l when T is small, the differ- 



ence between the SK model and the REM becomes trans- 
parent. The inverse participation ratio is close to unity 
in both cases and the main contribution comes from the 
first excited state for the SK model and from the ground 
state for the REM. We also see from the distribution 
function of I that / is always larger than 1/N which is 
consistent with results of the analysis of xi 10 ) f° r the SK 
model. Thus, the gap analysis above is justified for the 
SK model. On the other hand, the contributions of the 
excited states are very small and the classical analysis by 
static approximation is justified for the REM. 

We find that the linear, spin-glass, and nonlinear sus- 
ceptibilities diverge at different points, namely, T\, T sg , 
and r„i respectively. r sg corresponds to the phase tran- 
sition point. Although our system size is not sufficiently 
large, we roughly estimate that the extrapolation to the 
N — » oo limit gives r sg /J ~ 1.4, which is consistent with 



the previous numerical results llj, ll2|. In our numeri 



cal study, T n \/J ~ 1.6 is somewhat larger than r sg , and 
T\/J ~ 0.9 is close to the point where the average energy 
gap is zero. These behaviors can be verified by study- 
ing the distribution function of the susceptibility P x (x)- 
The power-law behavior P x (x) ~ X s with s < 1 at 
large x gives the divergence of the average susceptibility. 
This distribution is roughly the same as the power-law 
distribution of the gap since xi ~ for example. Our 
detailed study will be reported elsewhere. 

Although our result is very similar to that obtained in 
finite-dimensional systems, we cannot attribute the real- 
ization of the small gap to the presence of ordered clus- 
ters. This is because no notion of the spatial dimension 
exists for the present infinite-ranged model. As we men- 
tioned earlier, the fluctuation in extra imaginary time di- 
rection can induce the quantum Griffiths-like singularity 
where an ordered cluster is randomly formed in its direc- 
tion, which is similar to that obtained in low-dimensional 
classical systems. As in the finite-dimensional case [l|, 
by taking into account a rare event, we find a power-law 
behavior of the gap distribution function as 



P(A) 



r 5 A- A 



oe 



A fc , (7) 



where c, k, and Ao are constants and r is the imag- 
inary time. Instead of spatial directions in the finite- 
dimensional case [l|, a small gap is caused by a large 
imaginary time r with an exponentially small probabil- 
ity. From the facts that the randomness is absent in the 
imaginary time direction and Eq. ([7} is relevant only for 
the low-temperature limit of quantum systems, it follows 
that the effect should be distinguished from the Griffiths 
singularity. Since the present effect is present not only 
for the infinite-dimensional system but also for finite- 
dimensional ones, it is interesting to closely study how 
both effects are distinguished in terms of physical quan- 



tities in finite-dimensional systems. 

In conclusion, we have discussed quantum spin-glass 
transitions in terms of the energy gap. In the SK model, 
our analysis implies that three susceptibilities behave dif- 
ferently from each other and diverge at different points, 
which must be confirmed by further studies. The result 
can be characterized by a power-law behavior of the dis- 
tribution, which we attribute to the quantum Griffiths- 
like effect. In contrast to the classical Griffiths singu- 
larity, rare configurations of the quantum fluctuations 
play a significant role in the phase transition. On the 
other hand, the same mechanism cannot be applied to 
the REM. We have discussed conditions for this mecha- 
nism to be applied by studying several quantities. It is 
important to state the conditions in a more general way, 
which will be part of our future investigation. 
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